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Abstract 

The transcription factor nuclear factor kappa-B (NFkB) is a key regulator of pro-inflammatory and pro-proliferative processes. 
Accordingly, uncontrolled NFkB activity may contribute to the development of severe diseases when the regulatory system 
is impaired. Since NFkB can be triggered by a huge variety of inflammatory, pro-and anti-apoptotic stimuli, its activation 
underlies a complex and tightly regulated signaling network that also includes multi-layered negative feedback 
mechanisms. Detailed understanding of this complex signaling network is mandatory to identify sensitive parameters that 
may serve as targets for therapeutic interventions. While many details about canonical and non-canonical NFkB activation 
have been investigated, less is known about cellular IkBcx pools that may tune the cellular NFkB levels. IkBcx has so far 
exclusively been described to exist in two different forms within the cell: stably bound to NFkB or, very transiently, as 
unbound protein. We created a detailed mathematical model to quantitatively capture and analyze the time-resolved 
network behavior. By iterative refinement with numerous biological experiments, we yielded a highly identifiable model 
with superior predictive power which led to the hypothesis of an NFKB-lacking IkBcx complex that contains stabilizing IKK 
subunits. We provide evidence that other but canonical pathways exist that may affect the cellular IkB« status. This 
additional kBoclKKy complex revealed may serve as storage for the inhibitor to antagonize undesired NFkB activation under 
physiological and pathophysiological conditions. 
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Introduction 

The nuclear transcription factor kB (NFkB) family consists of 
five DNA-binding proteins (p65, p50, p52, cRel, RelB) that 
differentially modulate gene transcription. Since NFkB activation 
is involved in many cellular processes including inflammation, 
proliferation, angiogenesis and anti-apoptosis, only transient 
expression of the responsive genes ensures proper function of 
living cells [1]. Impairment of the regulatory system may 
contribute to malignant transformation, invasion and metastasis 
[2]. Consequently, NFkB controls its own activity by initiating 
negative feedback mechanisms including transcriptional up- 
regulation of cellular inhibitors [3,4]. Two independent pathways 
have been described to induce NFkB activation. While the non- 
canonical pathway determines the activity of p50:p52 heterodi- 
mers, the more prominent canonical pathway controls activation 
of p65:p50 subunits. In un-stimulated cells NFkB (p65:p50) resides 
inactively within the cytosol, bound to its inhibitor IkBoc, which 
covers its nuclear localization signal [5] . Following canonical signal 
transduction, NFkB activation is triggered via the IkB kinase 
complex (IKK), which consists of two catalytic subunits, IKKa 
and IKK(3, as well as the regulatory subunit IKKy. Site specific 
phosphorylation of IKKfi mediates downstream phosphorylation 
of IkBoc at two Ser residues, serving as a signal for poly- 
ubiquitination and proteasomal degradation of the inhibitor. 



Liberated NFkB subsequently translocates into the nucleus to 
serve its function as a transcription factor which also - and most 
importandy - includes induction of negative feedback regulation 
via IkBoi re-synthesis [6] . The network of interaction, however, 
leading to activation, inhibition or post-activational attenuation of 
NFkB is very complex and can be influenced by changes in signal 
transduction as well as through specific modifications of the 
molecules involved. Since dysregulation of NFkB plays a major 
role in the development of various diseases, a number of 
mathematical models have been implemented to analyze the 
non-linear dynamical behavior of this regulatory network [7,8]. 
These models contributed to the analysis of the role of negative 
feedback loops [9-1 1], the description of the overall input-output 
behavior of the pathways [12-14], the understanding of the 
integration from multiple input signals [14,15], as well as to the 
identification of sensitive parameters [16]. However, the pattern of 
NFKB-induced gene expression may significantly change depend- 
ing on the cell type, the intracellular protein-protein interaction 
and the physiological context. Following this line, previous studies 
revealed canonical NFkB responses to dramatically change in cells 
being exposed to DNA-damaging agents, including ultraviolet-B 
(UVB) radiation. In particular, Interleukin- 1 (IL-1) stimulation was 
shown to protect epithelial cells from death ligand-induced 
apoptosis via NFKB-dependent up-regulation of anti-apoptotic 
genes. When co-irradiated with UVB, instead, IL-1 driven and 
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Author Summary 

In unstimulated cells, the transcription factor NFkB resides 
in the cytosol bound to its inhibitor IkBoc. Canonical 
activation of NFkB by numerous stimuli leads to protea- 
somal depletion of IkBoc, thereby liberating NFkB to 
translocate into the nucleus to induce transcription of 
genes leading to proliferation, angiogenesis, metastasis, or 
chronic inflammation. Consequently, only transient activity 
needs to be warranted by immediate NFicB-dependent 
induction of negative regulatory mechanisms, including 
up-regulation of its inhibitor IkBoc. Resynthesized IkBoc 
consequently terminates NFkB activity by binding to its 
nuclear localization sequence. However, under physiolog- 
ical or pathophysiological conditions, random NFkB 
activation may occur, which needs to be avoided in order 
to guarantee proper cellular function. Using detailed 
dynamical modeling, we have now identified an additional 
IkBoc containing complex to exist in un-stimulated cells 
which lacks NFkB but includes IKKy (kBoclKKy complex). 
This additional IkBoc is not depleted from cells in the 
canonical fashion and may therefore serve as a cellular 
backup to avoid random NFkB activation. 



NFKB-dependent repression of anti-apoptotic genes caused 
enhancement of UVB-induced apoptosis [17]. As a prerequisite, 
nuclear persistence of NFkB was shown to be facilitated via UVB- 
induced inactivation of the catalytical subunit of Ser/Thr 
phosphatase PP2A, causing chronic IKKP activation and subse- 
quent phosphorylation and proteasomal degradation of resynthe- 
sized IkBoc. Both modifications in concert were shown to drive the 
pro-apoptotic properties of a classical non-apoptotic protein 
[3,18]. Using a systems biological approach we could show that 
not only PP2A inactivation but also global translational inhibition 
appeared to be involved in preventing IkBoc recurrence. Above 
this, translational inhibition was shown to induce IkBoc depletion 
in cells irradiated with UVB alone, indicating both mechanisms to 
individually mediate UVB-dependent responses [15,19]. 

Similar to UVB, co-stimulation of cells with IL-1 and the 
tyrosine phosphatase inhibitor orthovanadate (OVA) induced 
tyrosine kinase cSrc-mediated inactivation of PP2A. Again here, 
chronic phosphorylation of IKKP and consequently inhibition of 
IkBoc recurrence provided for sustained NFkB activation [18]. In 
the present study we aimed to understand which additional 
regulatory mechanisms may exist to prevent unwanted NFkB 
activation under physiological conditions. We therefore analyzed 
OVA-dependent IkBoc depletion as a tool to identify additional 
IkBoc sources within the cell by iterative model refinement in 
combination with model inspired experimentation. The extended 
model predicted alternative non-canonical iKBoc-degradation to 
occur without affecting NFkB activity. Respective experimental 
design finally revealed a yet unknown lKBoc:IKKy complex to 
exist, which might serve as a backup for negative feedback 
regulation of NFkB. 

Results 

OVA treatment alone causes incomplete IkBoc depletion 

Stimulation of cells from the epithelial cell line KB with IL- 1 
caused canonical degradation of the NFkB inhibitor IkBoc via 
phosphorylation of the upstream kinase IKKP at Ser 1 77/181. 
Perfectly in line with the phosphorylation pattern of IKKp 1 , IkBoc 
was completely degraded 30 min after IL-1 stimulation and 
became resynthesized after 2 h when phosphorylation of IKKP 



started to descent (Fig. 1A). Co-treatment of cells with IL-l+OVA 
instead accelerated phosphorylation of IKKP causing early 
phosphorylation and degradation of IkBoc (Fig. IB, 15 min). 
Strikingly, re-accumulation of IkBoc was completely inhibited 
under these conditions even though the phosphorylation pattern of 
IKKP at later time points (2 h; 4 h) remained largely unchanged 
compared to IL-1 treated cells (Fig. 1A). This strongly indicated 
that other OVA-driven mechanisms may superimpose IL-1- 
mediated canonical IkBoc degradation at later time points. In 
accordance with this assumption we could detect partial, almost 
linear IkBoc depletion at 4-8 h after OVA only treatment, 
following a much slower kinetics than canonical IL- 1 -driven IkBoc 
degradation (Fig. 1C). Since IKKP phosphorylation did not seem 
to play a major role in delayed OVA-dependent IkBoc depletion 
we next examined whether transcriptional or translational 
inhibition might be involved. 

IkBoc depletion is independent of transcriptional and/or 
translational inhibition 

Performing RT-PCR analysis we revealed the IkBoc mRNA 
level to remain completely unchanged, even up to 8 h after OVA 
treatment, while being up-regulated 1 h after canonical IL-1 
treatment, as a positive control (Fig. 2A). Application of the 
transcription inhibitor actinomycin D (ActD) strongly implied 
transcriptional alterations not to be involved in OVA-induced 
IkBoc depletion. While OVA treatment alone induced a moderate 
and incomplete reduction of the IkBoc protein (Fig. 2B) but not the 
respective mRNA (Fig. 2C) over time, transcriptional inhibition by 
ActD caused pronounced inhibition of both the mRNA and the 
protein level of IkBoc, respectively. Of note, co-application of OVA 
and ActD further enhanced depletion of IkBoc protein without 
additively affecting the transcription level (compare Fig. 2B + 2C). 
Results indicated that OVA-induced IkBoc depletion is facilitated 
at the protein level, independent of transcriptional regulation. An 
analogous IkBoc protein pattern was obtained when inhibiting 
translation by addition of cycloheximide (CHX). While individual 
treatment with either CHX or OVA caused IkBoc reduction over 
time, co-application of both substances additively enhanced loss of 
IkBoc (Fig. 2D). These data strongly support the concept that 
OVA-mediated IkBoc depletion is independent of translational 
inhibition but might be caused by activation of upstream signaling 
pathways apart from canonical NFkB signal transduction. 

OVA-induced IkBoc depletion neither follows the 
canonical pattern nor causes NFkB activation 

Canonical NFkB activation is well known to involve Serl77/ 
181 phosphorylation of IKKP, followed by Ser32/36 phosphor- 
ylation of IkBoc as a prerequisite for its proteasomal degradation. 
Stepwise documentation of canonical NFkB activation by 
Western-blot analysis and electro mobility shift assay (EMSA) 
revealed that OVA-induced IkBoc depletion does not follow the 
canonical pattern. While IL-1 -induced total IkBoc degradation 
occurred as a fast process being completed after 15 min, OVA- 
induced subtotal IkBoc depletion followed a much slower kinetics, 
and did not involve classical IKKP phosphorylation (Fig. 3A). 
Illustrating canonical IkBoc degradation by addition of the 
proteasome inhibitor MG132, which stabilizes phosphorylated 
IkBoc, revealed only canonical IL-1 stimulation to cause an IkBoc 
shift, while OVA treatment did not, but still caused IkBoc depletion 
over time (Fig. 3B). Correspondingly, OVA depleted not only 
IkBol-wI but also the Ser32/36Ala mutant which lacks the IKKP- 
dependent phosphorylation sites and can therefore not be 
degraded in the canonical fashion (Fig. 3C). Interestingly, EMSA 



PLOS Computational Biology | www.ploscompbiol.org 



2 



March 2014 | Volume 10 | Issue 3 | e1 003528 



Identification of New IkBc< Complexes 





kBa 



a-tubulin 



Figure 1. OVA promotes delayed depletion of IkBu independent of IKKp phosphorylation. A: KB cells were left untreated or pre-treated 
with OVA (1 mM) for 1 h followed by stimulation with IL-1 (10 ng/ml). At the indicated time points following IL-1 stimulation the cellular status of P- 
IKKp (Ser1 77/1 81), IkBoi and IKK(3 was determined by Western-blot analysis. Equal loading was monitored with an anti-oc-actin antibody. B: Cells were 
left untreated, stimulated with IL-1 or OVA for the indicated time points and the cellular IkBoi status determined by Western-blot analysis, a-tubulin 
served as a loading control, ***p<0.001. 
doi:10.1371/journal.pcbi.1003528.g001 



revealed no significant nuclear translocation of NFkB to occur 
upon delayed OVA-induced IkBoi depletion (Fig. 3D), indicating 
other than NFtcB-bound IkBoi pools to exist within un-stimulated 
cells. In order to identify this additional IkBoc pool we integrated 
all experimental data into detailed dynamical modeling using 
ordinary differential equations. 



kBa co-exists within different protein complexes 

Based on our previous NFkB signaling model [15] we generated 
an extended mathematical model (variant M-l, Fig. 4A). Since 
OVA-mediated IkBoi depletion appeared to be independent of 
both, transcriptional/translational inhibition and proteasomal 
degradation, we assumed proteases to be involved in this process. 
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Figure 2. OVA-dependent IkBh depletion is mediated at the post-translational but not at the transcriptional level. A: Cells were left 
untreated, stimulated with IL-1 (10 ng/ml) or OVA (1 mM). At the indicated time points transcriptional regulation of IkBoc was monitored by RT-PCR, 
with GAPDH serving as a loading control. B: Cells were pre-treated or not with ActD (5 |ig/ml) for 1 h and subsequently stimulated or not with OVA 
(1 mM) as indicated. Protein level of IkBoc was determined by Western-blot analysis with a-tubulin as loading control, and C: transcription level of 
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IkBoc by RT-PCR with GAPDH as loading control. D: Cells were pre-treated or not with CHX (5 ug/ml) for 1 h and subsequently stimulated or not with 
OVA (1 mM) as indicated. Protein level of kBa was determined by Western-blot analysis with a-tubulin serving as a loading control, **p<0.005; *p< 
0.05. 

doi:1 0.1 371 /journal.pcbi.1 003528.g002 



It is known that free as well as NFKB-bound IkBoc can be degraded 
by proteases in a proteasome- as well as in a lysosome-independent 
manner [20,21]. To monitor the role of OVA within the entire 
NFkB signaling network, we also included OVA-dependent Src 
phosphorylation - which subsequently leads to inactivation of 
PP2A and chronic IKK(5 activation - into the model. This part of 
the signaling pathway, however, exclusively plays a role at early 
time points during IL-1 dependent canonical IkBoc degradation 
[18]. Successful parameter fitting was performed based on a huge 
variety of experimental data comprising multiple stimulations: 
OVA-, ActD — , ActD+OVA, CHX+OVA-stimulation (Fig. 1-3), 
IL-1 + OVA-stimulation [18], as well as data describing IkBoc 



depletion upon IL-1-, UVB-, IL-1+UVB-, CHX—, IL-1+ 
UVB+MG132-stimulation collected before [15,19]. All experi- 
mental data could be fitted well (Fig. SI). In particular, the good 
reproducibility of OVA experiments indicated that addition of an 
iKBoc-degrading protease is sufficient to reproduce OVA-induced 
IkBoc depletion without activating NFkB. The model thereby 
predicted a pool of free IkBoc which is eliminated by a putatively 
OVA-activated protease, but leaves NFKB-bound IkBoc almost 
unaffected. In the best fit scenario of knock down simulation 
studies depletion of free IkBoc showed major relevance while the 
depletion of NFkB bound IkBoc is negligible (Fig. 4B). Of note, 
OVA-induced PP2A inactivation is only important in IL-1 
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Figure 3. OVA-induced IkBu depletion does not follow the canonical NFkB activation pattern. A: Cells were stimulated with IL-1 (10 ng/ 
ml) or OVA (1 mM) for the indicated time points. Status of p-IKKp and kBa was determined by Western-blot analysis. B: Cells were pretreated with 
the proteasome inhibitor MG1 32 for 1 5 min. Subsequently cells were stimulated as in A and the IkBoc status determined by Western-blot analysis. C: 
KB cells stably expressing the kBa super-repressor variant Ser32/36Ala were stimulated as cells in A, and the kBa status determined by Western-blot 
analysis. D: Cells were stimulated as in A, nuclear protein extracts were isolated and subjected to EMSA using an NFKB-specific consensus sequence. 
doi:1 0.1 371 /journal.pcbi.1 003528.g003 
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Figure 4. IkBcj co-exists in different cellular complexes. A: Schematic representation of model variant M-1 . Input (IL-1, UVB, OVA, MG132 (MG), 
ActD, CHX) is given in orange and red, respectively. Reactions and variables adopted from our previous model (Witt et al, [15]) are shown in grey. B: 
The red line represents simulation data of the best fit (reference). The impact of protease-mediated degradation of free IkBoc is depicted in orange 
and determined by setting the rate constant of this reaction (a 4 ) to zero. Likewise, the rate constant of protease-mediated degradation of NFkB 
bound IkBoc (a 5 ) was set to zero (green line). Experimental data and standard deviation is shown in blue. C: KB cells were left untreated, pre- 
stimulated with zVAD (20 |jJVI), CATI-1 (50 |jJVI) or calpastatin (1 |jJVI), for 1 h followed by OVA treatment for 8 h. Protein level of IkBoc was determined 
by Western-blot analysis with oc-tubulin serving as a loading control. D: Cells were left untreated, stimulated with IL-1 for 15 min or OVA for 8 h. 
Subsequently protein extracts were subjected to immunoprecipitation with an antibody specific for the p65 subunit of NFkB. The amount of co- 
precipitated IkBoc was determined by Western-blot analysis compared to input protein levels. Bars represent IkBoc levels that are normalized to the 
respective p65 levels. E: Lysates from unstimulated cells were subjected to size exclusion chromatography. Individual fractions were analyzed via 
Western-blotting, using antibodies against IkBoc, p65 and p50. F: IKKoc and IKKB, respectively, were transiently knocked down using siRNA, with siLacZ 
serving as a negative control. Subsequently cells were stimulated with OVA for 8 h and the IkBoc status, as well as IKKoc and IKKP knock down 
determined by Western-blot analysis. oc-tubulin served as loading control. 
doi:10.1371/journal.pcbi.1003528.g004 



stimulated cells and is therefore also irrelevant for OVA only 
stimulation following a much slower kinetics. 

Still, the model predictions had to deal with two obstacles: 
Firstly, due to its immense instability only a very small pool of free 
IkBoc seems to exist within the cell at all [22]. Secondly, proteolytic 
cleavage of IkBoc at least by three major groups of proteases: 
caspases, calpains and cathepsins could be excluded by the use of 
specific inhibitors (Fig. 4C). To follow up the model based 
hypothesis of a stably existing IkBoc form that is not bound to 
NFkB, we immuno-precipitated the NFkB p65 subunit from 
whole cell lysates and checked the levels of co-precipitated IkBoc in 
un-stimulated versus OVA-treated cells. The amount of IkBoc 
remained largely unchanged, whereas IkBoi was absent in cells 
stimulated with IL-1, due to complete canonical degradation 
(Fig. 4D). These data strongly supported the assumption that only 
IkBoc which is not bound to NFkB is depleted in an OVA- 
dependent fashion. To investigate whether depleted IkBoc in fact is 
free or bound to other cellular components but NFkB, we 
conducted size exclusion chromatography. Indeed, IkBoc appeared 
to exist in at least three different forms in un-stimulated cells. 
While only a minimal fraction seemed to refer to unbound IkBoc 
eluting at a size of 44 kDa, surprisingly no complex exclusively 
consisting of IkBoc and NFkB (p65:p50) seemed to exist. Instead, 
complexes of higher molecular weight containing IkBoc, p65 and 
p50 (NFkB) showed up at sizes ranging between 300 kDa and 
490 kDa, indicating to incorporate other proteins as well. In 
addition, different aggregates ranging between 100 kDa and 
300 kDa in size appeared which did not contain any NFkB 
components, implying IkBoc to also form complexes with other 
proteins (Fig. 4E). To investigate whether IKKs might be involved 
in stabilizing IkBoc we knocked down the two catalytic subunits 
IKKoc and IKKP by RNA interference. While IKKoc knock down 
had no effect on IkBoc depletion, knock down of IKKP seemed to 
enhance loss of IkBoc, implying that binding to IKKP stabilizes 
IkBoc which is not bound to NFkB (Fig. 4F). The fact that model 
variant M-1, and a simple extension by adding a second 
(competing) IkBoc complex formation (variant M-2, Fig. S3), 
respectively, could not fully explain the IKK knock down data 
(Fig. S2 and S4) argues for an additional IkBoc-IKK interaction 
and calls for a more in depth analysis of the IkBoc complexes. 

Novel cellular IkBoc complexes containing IKK 
components but no NFkB 

To stress whether IKK components play a role in IkBoc complex 
formation we first generated model variant M-3 assuming that 
IKK is permanently bound to IkBoc which lacks NFkB (Fig. 5A). 
According to the chosen model structure, binding of IKK to IkBoc 
results in the formation of an IkBocIKK as well as an 
NFkB:IkBoc:IKK complex. As phosphorylated IKK (IKKp) 
initiates degradation of IkBoc, the IKKp containing complexes 



are highly unstable and dissociate into free IKKp and free NFkB 
which translocates into the nucleus. 

The simulation results of this model variant revealed a very 
good reproducibility of all experimental data including the OVA+ 
IKK knock down experiment (Fig. 5B). 

When exploring the model predicted steady state levels in the 
un-stimulated IKK knock down setting (Fig. S5), very low IKK 
yielded decreased formation of both IKK containing complexes 
IkBocIKK and NFkB:IkBoc:IKK. As a consequence elevated levels 
of free and NFKB-bound IkBoc (NFkB:IkBoc) were predicted. 
Consequendy, constitutive degradation of NFkB:IkBoc results in 
enhanced NFkB activation and further increased levels of free 
IkBoc (Fig. S5A). Starting from these changed steady state levels in 
the IKK knock down setting - especially taking the increased level 
of free IkBoc into account - the OVA-mediated degradation of 
IkBoc could be reproduced accurately (Fig. S6). 

In this final model OVA mediated the degradation of an NFkB- 
free IkBoc complex (IkBocIKK), thus preventing activation of 
NFkB and IkBoc mRNA synthesis (Fig. 5C and Fig. 5D). 

In summary, the mathematical model based analysis clearly 
proposes the existence of IKK containing IkBoc complexes lacking 
NFkB, supporting the hypothesis that IKKP stabilizes "unbound" 
IkBoc. Size exclusion chromatography clearly revealed that the 
high molecular weight complexes (300-490 kDa) in un-stimulated 
cells consisted of at least IkBoc, NFkB (p65:p50), and all three IKK 
components IKKoc, IKKP and IKKy. Most importantly, IKKy 
eluted together with IkBoc at 1 00 kDa, perfectly matching the size 
of a heterodimeric complex (±94 kDa). Of note, additional 
complexes containing exclusively IKKP and IKKy (±140 kDa) 
but not IKKoc appeared to be formed (Fig. 5E). 

Treatment with OVA for 8 h did not significandy change the 
composition of the high molecular weight complexes. In contrast, 
it caused pronounced depletion of IkBoc from the lKBoc:IKKy 
complex, while IKKy seems to randomly distribute over numerous 
fractions. This clearly indicates dissociation of this heterodimeric 
lKBoc:IKKy complex to precede IkBoc depletion (Fig. 5F). The 
existence of an NFKB-lacking lKBoc:IKKy complex as well as the 
specific depletion of IkBoc from this particular complex following 
OVA treatment was confirmed by co-immunoprecipitation. While 
in unstimulated cells reasonable amounts of IkBoc but no NFkB 
were shown to bind to IKKy, the level of IkBoc decreased 
significandy upon treatment of cells with OVA for 8 h. No IkBoc 
was found bound to IKKy in IL-1 stimulated cells, due to 
complete canonical degradation (Fig. 5G). This is perfectly in line 
with simulation results derived from the final model variant M-3 
suggesting that degradation of NFKB-free but IKK-bound IkBoc is 
responsible for the partial IkBoc depletion in response to OVA. 

To assess the functional relevance under physiological condi- 
tions, we investigated IkBoc depletion in response to stimuli being 
capable to induce both, NFkB activation and apoptotic cell death 
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Figure 5. NFkB free licBa complexes contain subunits of the IKK complex. A: Schematic representation of model M-3 (final model). State 
variables are depicted in blue, inputs (IL-1, UVB, OVA, MG132 (MG), ActD, CHX) in orange and red, respectively. Reactions and variables adopted from 
our previous model (Witt et al, [15]) or model variant M-1 or M-2 are shown in grey. B: Overall simulation results of the best fit. The 127 data points 
and standard deviations are depicted in blue, the red line represents simulated time course of the respective state variable. The overall % 2 value is 
1 00.6. C: Reference (red line) represents simulation data of the best fit. Setting the rate constant of protease-mediated degradation of free IkBoe (a 4 ) to 
zero reveals the impact of this reaction on observed IkBoc degradation (orange). The impact of protease-mediated degradation of IKK bound IkBoc (p,) 
is depicted in green. D: Simulated time course of the IkBoc mRNA level after 8 h of OVA treatment is depicted in red. Corresponding experimental 
data and standard deviation is shown in blue. Lysates from E: unstimulated or F: OVA treated cells were subjected to size exclusion chromatography. 
Individual fractions were analyzed by Western-blotting, using antibodies against IkBoi, p65 and p50, IKKoc, IKKF3 and IKKy. G: Cells were left untreated, 
stimulated with IL-1 for 15 min or OVA for 8 h. Subsequently IKKy was immunoprecipitated from protein extracts and the amount of co-precipitated 
NFkB (p65) and IkBoc determined by Western-blot analysis compared to input protein levels. H: KB cells stably expressing the IkBoi super-repressor 
variant Ser32/36Ala were stimulated as indicated. After 24 h cellular IkBoi status was monitored by Western-blot analysis with oi-tubulin serving as a 
loading control. 

doi:10.1371/journal.pcbi.1003528.g005 



[23,24]. Treatment with the death ligands TRAIL and FasL, 
respectively, exclusively caused canonical IkBoc degradation, being 
indicative by sparing the IkBoi super-repressor variant from 
depletion. In contrast, irradiation of cells with UVB additionally 
resulted in IkBoc depletion independent of the canonical pathway - 
represented by degradation of both, the endogenous and the super- 
repressor variant (Fig. 5H). These data strongly support the notion 
that other but canonical pathways exist that may affect the status 
of IkBoc within the cell. Thus, we have uncovered an additional 
IkBoc complex to exist in un-stimulated cells that might serve as 
storage to antagonize random or undesired NFkB activation and 
consequently ensures proper cellular function. 

Discussion 

In the past decades activation of NFkB has exclusively been 
attributed to either canonical or non-canonical signal transduction 
pathways. Both pathways are tightly regulated by a series of 
phosphorylation and ubiquitination events that cause nuclear 
translocation of distinct NFkB family members. Canonical signal 
transduction basically accumulates at two well described protein 
complexes, namely IKK and NFkB:IkBoc. Accordingly, stimula- 
tion of cells with the pro-inflammatory cytokine IL-1 causes 
downstream activation of the IKK complex, in particular the 
catalytic subunit IKKP, to mark IkBoc for proteasomal degrada- 
tion. Released NFkB in turn triggers resynthesis of its inhibitor in a 
negative regulatory feedback loop. Most recently, a number of 
ubiquitin ligases including TRAF molecules and the lubac 
complex have additionally been implemented in canonical NFkB 
activation [25,26], while de-ubiquitinases like A20 and CYLD 
serve a well-known function in negative feedback regulation 
[27,28]. In the present study we provide evidence that besides the 
well-known components of canonical NFkB signaling, alternative 
IkBoc containing complexes exist within the cell that might 
indirectly contribute to NFkB regulation. Short term inhibition of 
IkBoc resynthesis following IL-l+UVB and IL-l+OVA treatment, 
respectively, is due to abrogation of canonical negative feedback 
regulation via inhibition of PP2A leading to chronic IKKP 
activation [3,19]. At later time points, however, canonical 
feedback regulation seems to be superimposed by an IKK(3- 
independent OVA- driven mechanism that follows a slower 
kinetics and only causes partial IkBoi depletion. A similar 
observation could previously be made in UVB-irradiated cells, 
showing slow and subtotal IkBoi depletion resulting in only 
moderate and delayed NFkB activation [3]. In the present study 
OVA-induced delayed IkBoc depletion appeared to follow a 
pattern different from canonical NFkB activation, because super- 
repressor variants of IkBoi could be depleted as well. This implies 
IkBoi depletion which does not follow the canonical pattern to 
serve an important, yet unknown function. Although a number of 
alternative ways to proteolytically cleave IkBoc have been 



described in the literature [20,21,29,30] inhibition of the major 
cellular protease families, could be ruled out. More recently an 
alternative proteasome independent mechanism called PIR has 
been described to enhance IkBoc turnover in B-cells, however, PIR 
resulted in constitutive p50:cRel activation in those cells and may 
therefore play a different role [21,31]. Still, integration of an 
OVA-activated protease into all our mathematical models that is 
able to deplete IkBoc- different from the canonical mechanism - 
was shown to nicely reproduce the OVA induced IkBoc 
degradation with slow activation kinetics determined by a small 
rate constant (M-3: kp rot = 3.76e -7 s _1 , Table SI). Accordingly, 
significant degradation of IkBoc occurs at later time points (4 h- 
8 h) and may involve PIR-like mechanisms. 

Reproducing OVA-induced IkBoc depletion without NFkB 
activation in model variant M-1 predicted high levels, 31%, of free 
IkBoc (0.041 |xM out of 0.135 u,M) to exist within the cell (Table 
S2), whereas only 10% to 15% of free IkBoc is supposed to exist 
[32,33]. Remarkably, inclusion of putative NFKB-lacking IkBoc 
complexes nicely scaled down the level of free IkBoc to 13% (M-2) 
and 8% (M-3) respectively which now is in very good accordance 
to the literature values, matches mathematical models of other 
groups [11,22,34], and could also be verified by size exclusion 
chromatography. 

According to previous studies that revealed IKKa and IKKP to 
form high molecular complexes that contain IkBoc as well as NFkB 
components [35,36] model variant M-3 predicted IkBoc to form 
both, stable IkBocIKK as well as NFkB:IkBoc:IKK complexes. 
With this final model the entire data of numerous experiments 
could be reproduced with a good fit quality. Finally formation of 
an iKBocIKKy complex could be verified experimentally by gel 
filtration analysis as well as co-immunoprecipitation. Expanding 
the scope of our previous model [15] by iterative model refinement 
revealed new insights into NFkB regulation and allowed to analyze 
the effect of an in silico knock down of the IKK complex on the 
OVA-mediated IkBoc depletion. The model predicts that knocking 
down IKKP enhances the level of an IkBoc compound that can be 
degraded by the proposed protease. In vitro neither IkBoc from the 
high molecular complex nor IKKP:IKKy bound IkBoc seems to be 
depleted after OVA treatment. Thus knocking down IKKP could 
result in decreased levels of NFkB:IkBoc:IKK as well as 
IkBoc:IKKP:IKKy and an elevated concentration of the Ik- 
Boc:IKKy complex. Due to larger amounts of IkBocIKK/, OVA- 
mediated degradation of this IkBoc component would explain why 
IKKP knock down enhances OVA-induced IkBoc depletion. 
Additionally, a constant basic IKKP-mediated IkBoc turnover 
exists in unstimulated cells [19], yielding in continuous NFkB- 
dependent resynthesis of IkBoc. In case of IKKP knock down, this 
process is completely abrogated consequendy causing subtle loss of 
IkBoc over time exclusively dependent on the individual half life in 
treated versus untreated cells (see also Fig. 4F). 
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Combining experimental methods and detailed dynamical 
modeling we could provide strong evidence for the existence of 
an NFKB-free, IKKy containing IkBoi complex presumably acting 
as a cellular backup pool to capture randomly released NFkB. 
Above this, our model introduces an IkBoc degradation pathway 
that is independent from canonical processes but is likely to 
influence the cellular status of IkBcz. Our final mathematical 
model (M-3) created here is able to reproduce a huge number of 
datasets comprising various intracellular proteins and a wide range 
of stimulation experiments to extensively reflect on the whole 
NFkB signaling network above individual top-down signaling 
pathways. Thus, our model can be used for further modeling 
approaches regarding NFkB regulation and may provide predic- 
tive potential for sensitive parameters that may serve as 
therapeutic targets in the future. 

Materials and Methods 

Cells and reagents 

The human epithelial carcinoma cell line KB (ATCC) was 
cultured in RPMI 1640, 10% FCS. Recombinant human IL-ip 
(R&D Systems, Wiesbaden, Germany) was applied at 10 ng/ml 
and Na-Orthovanadate (Sigma, Munich, Germany) at 1 mM. 
Actinomycin D and cycloheximide (Sigma) were added to cells at 
5 ng/ml, respectively. Specific protease inhibitors (Calbiochem, 
Darmstadt, Germany) were applied at 50 uM for the cathepsin 
inhibitor CATI-1, 1 uM for calpastatin and 20 uM for the pan 
caspase inhibitor zVAD Proteasomal inhibition was achieved by 
addition of 25 uM MG132 (Calbiochem). For Fas/CD95 receptor 
activation 0.5 ug/ml of an agonistic antibody (Immunotech, 
Monrovia, CA, USA) was used. Recombinant human iz-TRAIL 
protein, N-terminally fused to an isoleucine-zipper motif in order 
to constitutively build the trimerized active form [37] was kindly 
provided by Dr. Henning Walczak, Centre for Cell Death, Cancer 
and Inflammation, UCL, London and added at 100 ng/ml. UVB 
irradiation (300 J/m 2 ) was performed with TL12 fluorescent bulbs 
(290-320 nm, Philips). 

Immunoprecipitation and WB analysis 

Cells were lysed in lysis buffer (50 mM Hepes, pH 7.5; 150 mM 
NaCl; 10% glycerol; 1% Triton-X- 1 00; 1.5 mM MgCl 2 ; 1 mM 
EGTA; 100 mM NaF; 10 mM pyrophosphate, 0.01% NaN :i and 
Complete protease inhibitor cocktail; Roche, Mannheim, Ger- 
many) for 20 min on ice. Endogenous NFkB (p65) or IKKy were 
immune-precipitated using specific antibodies (sc-372; sc-8330 
Santa Cruz, Heidelberg, Germany) and A/ G-plus agarose (Santa 
Cruz) over night. Precipitates were analyzed by Western-blotting 
using antibodies against NFkB (F6, sc-8008, Santa Cruz), IkBoi 
(L35A5, CeU Signalling, Beverley, MA, USA) and IKKy (IMG- 
324A, Imgenex, San Diego, CA, USA). For WB analysis cells were 
lysed by addition of hot (95 °C) Laemmli buffer. 80 Ug protein 
extracts were subjected to SDS-PAGE and Western-blot analyses 
using antibodies against IkBoi, P-IKK(3-Serl77/181, IKK|3, p50 
(L35A5, 16A6, 2C8, #3035, Cell Signaling), p65 (sc-8008, Santa 
Cruz), IKKoc (556532, BD Biosciences), IKKy (IMG-324A 
Imgenex, San Diego, CA, USA), and oc-tubulin (DM1 A, 
Neomarkers, Fremont, CA, USA), using West-Pico or West-Dura 
(Pierce, Thermo Scientific, Rockford, IL, USA) chemiluminescent 
substrates. 

Semiquantitative RT-PCR analysis 

Total RNA was extracted from cells using GIT-buffer (4 M 
guanidinthiocyanate, pH 4.8; 0.3 M NaOAc; 1% N-lauroylsarco- 
sine; 0.2% P-mercaptoethanol) followed by phenol/chloroform 



extraction utilizing Phase Lock Heavy tubes (Eppendorf AG, 
Hamburg, Germany). Six ug of total RNA was reverse transcribed 
with an AMV Reverse Transcriptase kit (Promega, Mannheim, 
Germany). The following primers were used in a 20 ul reaction 
utilizing the RedTaq polymerase system (Sigma): 
GAPDH: 

F: 5'-GCCTCCTGCACCACCAACTGC-3'; R: 5'-CCCT- 
CCGACGCCTGCTTCAC-3' 
IkBcz: 

F: 5 ' -ACAGGAATTACAGGGTGC AGG-3 ' ; R: 5'-GAGA- 
AACTCCCTGCGATGAG-3' 

Plasmids and transfection of cells 

For ectopic expression of IkBoi -S32/36A 6, 5xl0 6 cells were 
transfected with 25 ug of the respective pcDNA3.1 -based 
construct by electroporation at 1200 uF and 250 V (EasyjecT- 
plus, Peqlab, Erlangen, Germany) in ice cold RPMI medium w/ o 
FCS. Transfection efficacy ranged between 70 and 80%. 

RNA interference 

1x10 cells were transfected with 0.5 pmol/ul siRNA knocking 
down IKKoe: GCAGAAGAUUAUUGAUCUATT or IKKfj: 
UUCAGAGCUUCGAGAAGAATT (Eurofins MWG, Ebers- 
berg, Germany) using Lipofectamin 2000 (Life Technologies, 
Darmstadt, Germany) according to the protocol. Cells were 
analyzed after 72 h. 

Electro mobility shift assay (EMSA) 

Following stimulation cells were harvested and nuclear proteins 
extracted as described before [38]. The NFkB consensus oligo 
nucleotide (sc-2505; Santa Cruz) was end-labeled using [y' 52 P] 
ATP and T4 polynucleotide kinase (MBI Fermentas, Ontario, 
Canada), followed by column-purification (QIAquick Nucleotide 
Removal Kit, Qiagen, Hilden, Germany). Binding reactions were 
carried out in a 20 ul volume containing 8 Ug nuclear protein 
extract in 5x binding buffer (20 mM HEPES, pH 7.5; 50 mM 
KC1; 2.5 mM MgCl 2 ; 20% (w/v) ficoll; 1 mM DTT), containing 
1 ug polyfdldC]; 2 Ug BSA, and 70.000 cpm of 32 P-labeled NFkB 
consensus oligo nucleotide for 20 min at RT. Samples were 
separated on a 4% native PAGE at 150 V for 2.5 h and detected 
by autoradiography. 

Size exclusion chromatography (FPLC) 

Cells were harvested in PBS+0.01% sodium acide and disrupted 
by sonication. 10 mg protein extract was fractionated on an 
agarose bead column (GE, Healthcare, Frankfurt, Germany). 
Fractionated proteins were precipitated in 50% TCA, washed 
twice with acetone and applied to SDS-PAGE and subsequent 
Western-Blot analysis. Proteins of known size (thyroglobulin: 
670 kDa, y-globulin: 158 kDa, ovalbumin: 44 kDa, myoglobin: 
17 kDa, cobalamin: 1.3 kDa) were used as a standard. 

Data processing 

Western-Blots analyses are presented as mean ± SD of 3 
independently performed experiments. Blots were evaluated 
densitometrically and the results normalized to the maximal 
measured value. A minimal relative standard deviation of 15% was 
always assumed. For statistical analysis student's t-test was 
performed. 

Mathematical modelling 

The mathematical models created in this study are based on our 
previous ODE model [15] comprising the following components: 
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IL- 1 -receptor (IL-1R), IL-1 -receptor ligand complex (ILRc), 
free IkBoc (IkBoc), NFkB bound IkBoc (NFkB:IkBoc), nuclear IkBoc 
(IkBocJ, nuclear NFkB (NFkBJ, IkBoc mRNA (IkBoc,), IKKP 
(IKK), phosphorylated IKKP (IKKp), Protein phosphatase 2A 
(PP2A). 

In our recent models IKK (and IKKp) is perceived as the IKK 
complex consisting of IKKot, IKK(3 and IKKy instead of IKKP 
alone. 

First of all we extended the basic model [15] by our previously 
proposed mechanism of OVA prolonging the IL-1 mediated 
activity of NFkB [18]. OVA treatment causes phosphorylation of 
the tyrosine kinase cSrc which in turn phosphorylates and thereby 
inactivates PP2A: 



j SRCp(t) = k om -(\ - SRCp(t))ova(t) 

d -PP2A(t)= -k„uv(t)-PP2A(t) 
at 

- k src -PP2A(tySRCp(t) + k ppdp -PP2Ap(t) 
~PP2Ap(t)=k m PP2A{tySRCp(t)-k ppdp -PP2Ap{t) 



Terms adopted from the model of Witt et al. [15] are 
highlighted in bold font. PP2Ap represents phosphorylated 
PP2A. Without loss of generality the total amount of the kinase 
cSrc (SRC) can be assumed to be 1 (cf. an analogous reasoning for 
IKK in [19]). Considering mass conservation, (l-SRCp(t)) then 
represents the Src kinase in its inactive (un-phosphorylated) state. 
The stimulus OVA is represented by the step function ova(t) having 
a value of 0 (absent) or 1 (present). The initial concentration of 
PP2A is set to 1 whereas PP2Ap is set to zero. In order to reach 
steady state conditions we started stimulation from steady state 
established after 240 hours of simulation of the un-stimulated 
system. 

Based on our recent results we assumed for the observed OVA 
mediated IkBoc degradation the involvement of a protease prot that 
is activated by OVA. This assumption is included in all of our 
model variants by adding the following term: 

^-prot(t) = k prM \\ -prot(t))-ova{t) 



In this equation prot represents the active form of the assumed 
protease and (l-prot(t)) the inactive state considering mass 
conservation. The total amount of prot is assumed to be 1 without 
loss of generality. In model variant M- 1 this protease is able to 
degrade free and NFkB bound IkBoc 

d 

— lKBu(t)=-a r (\ -mg(t)yiKKp(tyiKBot(t) 
at 

+ C4o'(l — chx{t))-{\ — uv(t)-uvinh)-lKBa,(t) 
- c 5a (l - mg(t)yiKBa(t) - h a lKBa(f) 
+ e\alKBoL„(i) — a<\-prot(t)-lKBa{t) 



— NFkB/ lKBu{t) = -ay{\ -mg(t)yiKKp(tyNFKB/lKBa{t) 

+ (aj iKBaJtyNFKBM))- r - 
k v 

-c 6a -(\ -mg(tyyNFKB/lKBa(t) 
-a 5 -prot(tyNFKB/lKBoc(t) 

Degradation of NFkB bound IkBoc by the protease in model 
variant M-l causes translocation of NFkB into the nucleus 
followed by the initiation of IkBoc mRNA synthesis: 

d 

— NFKB„(t) = k r -ay{\ -mg{t)yiKKp(t)-NFKB/ Ii<Bx(t) 

—a\ -IkBoCi, (t)-NFKB„ ( t) 

+ k v -c 6a -(l -mg(t)yNFKB/ IxBait) 

+k v -ayprot(tyNFKB/lKBu.(t) 

d 

— lKB<Xt(t) = (ci a -NFKB n (t) +C, rm sc)i 1 ~ ActD(t)) - C3alKBct t (t) 



The factor k v is used to compensate the different volumes of 
cytosol and nucleus (see Witt et al. [15]). 

ActD(t) represents a step function of the inhibitor of transcription 
actinomycin D (ActD) with a value of 1 (present) or 0 (absent). 

RNA polymerase II and basal transcription factors are known to 
be essential for transcription of genes [39,40] whereas specific 
transcription factors function as enhancer or repressor. Thus we 
assumed in M- 1 and M-2 a constitutive transcription rate for IkBoc 
mRNA synthesis (c fa „ JC ) that is independent of NFkB and is also 
included in the NFkB model of Hoffmann et al. [9] . 

Model variant M-l was fitted to time courses of NFkB, IKKP - 
P, IkBoc mRNA and/or IkBoc respectively gained from the 
following stimulation experiments: 

• IL-1 (IkBoc, IKKP -P, NFkB), IL-l+UVB (IkBoc, IKKp -P, 
NFkB), IL1+UVB+MG132 (IkBoc), UVB (IkBoc), cyclohexi- 
mide (CFIX) (IkBoc), UVB (IkBoc) [15] 

• IL-l+OVA (IkBoi, IKKP -P, NFkB) [18] 

• ActD (IkBoc), ActD+OVA (IkBoc), CHX+OVA (IkBoc), OVA 
(IkBoc, IkBoc mRNA), OVA+IKK knock down (IkBoc) 

The in silico knock down of IKK is realized by adding a factor 
named silKK to the model variant and setting its value to 0.93 
which reduces the initial concentration of IKK to the experimen- 
tally determined 7% of IKKP concentration measured in 
untreated cells: 

IKK(Q) = (\-siIKK(t)) 

In all stimulation experiments distinct from the IKKP knock 
down experiment the value of silKK is set to zero. 

In model variant M-2 we included the proposed IkBoc complex 
lKBa:Comp, consisting of IkBoc and components Comp distinct from 
NFkB: 
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IkBo I Compif) = api -lKBo(t)-Comp{t) — apu IkBo / Comp(t) 

—p 2 -prot(t)-lKBo./Comp(t) 

-pi IkBo./ Comp(i) IKKp(f) (\ -mg(t)) 



— Comp(i) = — api -lKBo.{i) ■ Comp(t) + apu IkBo / Comp(t) 

+P2'prot(t)lKBoi/Comp(t) 

+ p i -IxBa/ Comp(t) IKKp(t) ( 1 - mg(t)) 



In this model variant, the activated protease prot additionally 
degrades IkBoc from the IkBoc complex, lKBoc:Comp. Since this 
complex is part of the measured cellular IkBoc concentration the 
overall IkBoc concentration in the model, IkBoc,,/,,, becomes: 



IkBol (t) 

IkBcx obs = iKBo(t) + NFkB/ IkBo + — + IkBo/ ' Comp(t) 



The final model, variant M-3, (Fig. 5A) was designed to 
investigate the possibility of IKK being part of the proposed 
lKBo:Comp complex. We therefore extended model variant M-2 by 
binding reactions of IkBoc and IKK or IKKp and removed the 
component lKBo:Comp. In contrast to M-l and M-2 we fitted the 
start concentration of IKK by adding the parameter IKKstart: 



IKK(0) = {\ -siIKK(t)ylKKstart 



In addition we assumed that IKK as well as phosphorylated 
IKK is able to bind IkBoc in its free and bound form. 



dt 



lKBa(t) = - ap 3 IKK(t) lKBa(t) - ap x -IKKp{t)-lKBo(t) 

+ ap n IkBo /IKKp(t) + ap 33 IkBo / IKK(t) 
+ (1 — chx{t))-{\ —uv(t)-uvinh)-cn a -lKBa t {t) 
— csa - (1 — mg(t)) -lKBa(t) — i\ a -lKBo(t) 
+ e\ a -lKBa„(t) — a 4 -prot(t)-lKBo(t) 



— IKK(t) = - k p ILRc{t)IKK(t) + k dp -PP2A(t)- 



IKKp(t) 



dt 



K m + IKKp{t) 
- k pcomt -IKK{t) - a Pi I K Ba{i) IKK(t) 
- ap A -IKK(t)-NFKB/lKBx(t) + ap 33 IxBa /IKK(t) 
+ap w NFKB/ IkBo./ IKK(t)+pYprot(t)lKBa/ IKK(t) 



+ k pcomt -IKK(t) - api -IKKp(t)-lKBo(t) 

+ap xv lKBoi/ IKKp(t)-ap 2 -IKKp(t)-NFKB/ lKBa(t) 

+ apirNFiiB / IkBo./ IKKpit) 

+ a 2 -lKBa/IKKp(t)-(l -mg{t)) 

+ ayNFKB/lKBa/IKKp(t)i\ -mg(t)) 



—lKBa/IKK(t) = apylKBa(t)IKK(t)-apiylKB«./IKK(i) 
dt 

- k p ILRc{t) lKBa. /IKK(t) 

I k a -PRIAM- lKB *l IKR P^ 
+ k dp FFlA{t) Km + lKM/IKKp{t) 

-p\ prot(t)-lKBa/IKK(t) 

- k pconst -lKBx/IKK{t) 



— IicBtx/IKKp(t) = ap x lKBa(t) IKKp(t) -ap u -IkBo./ 'IKKp(t) 

+ k p ILRc{t)lKBa/IKK(t) 
- a r IkBol/ IKKp{t)-{\ -mg(t)) 
IkBol/ IKKp{f) 



k dp -PP2A(t)- 



dt 



K m + IkBol/ IKKp(t) 
+ k pcomt -lKBo/IKK(t) 

NFkB/ IkBol/ IKK(t) = ap A -IKK(t)-NFKB/I K Bo(t) 

- ap 44 -NFKB/I K Bo/IKK(t) 
- k p -ILRc(t)- NFkB/ IkBo /IKK(t) 
NFkB/ IkBol /IKKp(t) 



+ k dp PP2A(t) 



K n , + NFkB/ IkBo./ IKKp(t) 
k pC onst'NFKB/ IkBol/ IKK(t) 



dt 



NFkB I 'IkBo/ 'IKKp(t) = ap 2 -IKKp(t)-NFKB/ iKBo(t) 

-ap 22 - NFkB/ IkBo/ IKKp(t) 
-ay NFkB/ IkBo./ IKKp(t){\-mg{t)) 
+ k p -ILRc(t)-NFKB/ IkBo / IKK(t) 
NFkB/ 'IkBo/ 'IKKp(t) 



-k dp -PP2A 



K m + NFkB /IkBo/ IKKp 
+ kpconst' ^FkB/ IkBo /IKK(t) 
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In contrast to variants M- 1 and M-2 a constitutive IkB« mRNA 
synthesis {c tmm ^j in M-3 is not necessarily required to reproduce all 
experimental data since fitting M-3 without c transc to experimental 
data results only in a slight decrease of the fit quality (overall 
% 2 = 100.6 compared to 98.6). Thus parameter c tramc is not 
included in our final model M-3. This is in line with the models 
of Lipniacki et al. and Ashall et al. assuming that only nuclear 
NFkB initiates ItcBa mRNA synthesis [11,41]. 

Likewise, we did not integrate a constitutive degradation of 
IkBch in the NFkB:IkBoc:IKK complex analogous to constitutive 
degradation of IkBoc in NFkB:IkB« (c 6n ) since inclusion of this 
reaction did not improve fit quality. 

The entire ODE system of each model variant is shown in Text 
SI, S2, S3. 

For parameter estimation as well as for solving and analyzing 
the ordinary differential equation system we used the MATLAB 
(Mathworks) toolbox PottersWheel [42]. For optimization the % 2 
value of the following objective function was minimized by using a 
trust region approach: 

The % -value depends on the estimated parameter values. 
Variable y represents the i-th measured value whereas f states the 
simulated state value at time point i and is dependent on the 
parameter values p. The factor a ; represents the standard deviation. 

Besides the newly introduced parameters the fitted parameters 
of the previous model were also included in the parameter 
estimation with the same parameter boundaries (Table SI). 
Scaling parameters were included to take the lack of absolute 
values of the experimental data into account. 

Fit sequences were applied in the estimation process: the starting 
value of each parameter in a sequence was thereby calculated as 

p* =pW 

The Factor p represents the fitted parameter value of the 
currently best fit of the fit sequence (PottersWheel F3 routine) and 
the variable s determines the strength of disturbance where 
e~N(0,n). Four subsequent runs were performed with 100 fits each 
using n = 4,1,0. 1,0.01, respectively. 

In addition we performed an identifiability analysis of the final 
model M-3, using the top 10% of 300 fits with random starting 
conditions. Many parameters are well identifiable with a low 
relative standard Deviation (Table SI). Furthermore for some of 
the less identifiable parameters, linear or non-linear correlations 
with other parameters exist which indicate that combinations of 
these parameters are identifiable (Table SI). 

All models are provided as PottersWheel model datasets in the 
supplemental data. 

Supporting Information 

Dataset SI PottersWheel model file of model variant 
M_l. 

(TXT) 

Dataset S2 PottersWheel model file of model variant 
M_2. 

(TXT) 

Dataset S3 PottersWheel model file of model variant 

M_3. 

(TXT) 



Figure SI Overall simulation data of model variant 

M-l. The 125 data points and standard deviations are depicted in 
blue. Simulation data is shown in red with an overall % 2 value of 
97. 
(PDF) 

Figure S2 Overall simulation data of model variant M-l 
including the IKK knock down experiment. Fitting M-l 
additionally to the IKK knock down experiment increases the % 2 
value to 129. The 127 data points and standard deviations are 
depicted in blue, simulation data is shown in red. 
(PDF) 

Figure S3 Schematic representation of model M-2. State 
variables are depicted in blue, inputs in orange and red, 
respectively. Reactions and variables adopted from previous 
models (M-l or Witt et al. [15]) are shown in grey. 
(PDF) 

Figure S4 Overall simulation data of model M-2 
including the IKK knock down experiment. The model 
was fitted to 127 data points and revealed an overall % 2 value of 
119. Experimental data and standard deviation is depicted in blue. 
The red line represents simulation data of the best fit. 
(PDF) 

Figure S5 Simulated time courses of model variant M-3 
from start of integration with or without IKK knock 
down. The impact of the IKK knock down on steady state 
concentration of free iKBot (A), NFkB:IkBoc (B) and nuclear NFkB 
(C) is depicted in green. The simulated time courses without any 
input are shown in orange. 
(PDF) 

Figure S6 Impact of a 4 , the rate constant of protease- 
mediated degradation of free iKBa, after IKK knock 
down and OVA treatment for 8 h. Setting a 4 to zero reveals 
the influence of this reaction on observed IkBc* concentration 
(green). Reference (red) represents the simulation data of the best 
fit whereas experimental data and standard deviation is depicted 
in blue. 
(PDF) 

Table SI Description and correlation analysis of pa- 
rameters of model variant M-3. For calculation of the 
standard deviation (SD) and correlation analysis the best 10% of 
300 fits are used. 
(PDF) 

Table S2 Portions of IkB(x components. 

(PDF) 

Text SI ODE system of model variant M-l. 

(PDF) 

Text S2 ODE system of model variant M-2. 

(PDF) 

Text S3 ODE system of model variant M-3. 

(PDF) 
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